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INTRODUCTION 


1 .1 Scope and Purpose of Report 

This progress report covers activities during the third 
quarter, November 25, 1975 through February 25, 1976, for 
LANDSAT Investigation #23790. This investigation is funded 
for 19 months to develop techniques in Texas state agencies 
for using LANDSAT data to inventory and monitor coastal re- 
sources and uses. The General Land Office (GLO) is the Texas 
agency coordinating this investigation. Other participating 
agencies are the Bureau of Economic Geology (BEG), Texas 
Water Development Board (TWDB), and Texas Parks and Wildlife 
Department (TPWD). 

1.2 Summary of Work Performed 

During the third reporting period most of the accomplish- 
ments have dealt with refining techniques for extracting infor- 
mation from LANDSAT data. Specific accomplishments for the 
third quarter have been 1) to modify part of the detection and 
mapping (DAM) package so that printer maps, scaled and registered 
to U.S.G.S. 7 1/2 minute topographic maps, can be generated for 
the ADP classification results, 2) to develop a biological verifi- 
cation procedure for documenting plant species that occur in marsh, 
3) to apply new procedures for correlating the image interpretation 
line boundary map with the scaled and registered ADP display, 4) 
to initiate examination of a new area, test site 2, using proce- 
dures developed in test site 3, and 5) to summarize the interim 
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status of data acquisition, data analysis, and information display 
capabilities for the design of the LANDSAT- based monitoring 
system. 


2.0 PROBLEMS 

2.1 Technical Problems 

No technical problems relating to image interpretation 
procedures were encountered during this quarter. The techniques 
outlined in Section 3.3 of the second quarterly report {Jones, et al . , 
1975) have been applied to data from test site 2 (Scene 1289 - 
16261, 8 May 73) without the need for modification. 

Problems encountered in the ADP classification of the LANDSAT 
data centered around the gathering of training statistics discussed 
in Section 3.1.2. In clustering a test site, several classes are 
produced which have large standard deviations. These classes are 
small (less than }% of pixels included in training area), scattered, 
and usually associated with shoreline features at the land/water 
transition zone. Breaking such classes into their composite sub- 
classes to achieve tighter standard deviations has been difficult, 
time consuming, and generally unsatisfactory. While some sub- 
classes have been isolated, others have eluded detection because 
of their small size, and hence, remain unclassified. 

Another problem involves the look-up table in the ELLTAB 
classifier. Programming constraints currently limit the combined 
table size to 12,000 words. Due to the spectral complexity of 
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the coastal areas, this limit is often exceeded. Solutions being 
explored involve either 1) program modifications to allow in- 
creasing the table size, or 2) separate classification of land and 
water features to reduce the size of the table needed in a given 
run. 

2.2 Staffing Changes 

Bureau of Economic Geology fBEG) . Mr. Robert Baumgardner 
has replaced Mr. Samual Shannon as assistant to Dr, Robert Finley 
at the Bureau of Economic Geology. Mr. Baumgardner has a B. A. 
degree in zoology from the University of Texas at Austin and is 
currently working toward an M. A. degree in geology at this 
same school. 

Texas Parks and Wildlife Department CTPVJD) . Mr. Larry Lodwick, 
a staff research botanist in the Parks Division of TPWD has been 
assigned part-time from the Austin office. Mr, Lodwick's duties 
include 1) assisting George Clements by designing an appropriate 
field approach to document plant composition within the marshes, and 
2) consulting on the definition and correlation of coastal vegetation 
with LANDSAT marsh classes. 

3.0 ACCOMPLISHMENTS 

3.1 Examination of Test Site 3 (San Antonio Bay Areal 
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3.1.1 Area Description of Test Site 3 

Test site 3 overlays parts of San Antonio and Espiritu 
Santo Bays (Figure 1), which contain over 16,500 acres of salt 
marsh located mostly along the bayward shore of Matagorda Island, 
including the marshes at Pass Cavallo near Port O'Connor. Char- 
acteristic salt and brackish marsh plant species include 
Spartina alterni flora , Sparti na spartinae , Distichlis spjcata , 
Monanthochloe littoral is , Batis maritima , and Borrichia frutescens . 

Approximately 12,800 acres of brackish to fresh marsh is 
found in the upper reaches of San Antonio Bay, above Hynes Bay 
and Mission Lake on the delta of the Guadalupe River. Character- 
istic plants in these marshes include Phragmites cominunnis , 

Spartina patens , Spartina spartinae , and Juncus spp. 

The climate of this portion of the Texas Gulf Coast is 
subhumid, warni temperate. Average annual rainfall is 34 inches 
and the temperature ranges from about 20“ F to 90“ F or more during 
the summer months. Winters are usually very mild, having an 
average of only seven days of freezing or below freezing weather 
per annum (Texas Parks and Wildlife Department, 1975). 

The area is relatively rural, having only two major industrial 
plants in the county (such as Union Carbide near Mission Lake on the 
Austwell quadrangle) and a few small communities characterized by 
widely spaced buildings and large, grass-covered lots. These 
small communities, such as Austwell (population 284) and Port 
O'Connor (population 350) are devoted to farming, or sport and 
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ccniniercia’I fishing (The Dallas Morning News, 1972). The acreage 
of rangeland devoted to beef production is over twice that devoted 
to croplands in the area. The two major crops are rice and sorghum. 
Minor crops include cotton, soybeans, flax and sunflowers. The 
majority of the rangeland is unimproved, consisting of mixed native 
grasses, shrubs, and trees (Gilbert Hei deman, County Agricultural 
Extension Agent, personal communication, 1976). 

The areas of emphasis within test site 3 are those with the 
most extensive wetlands, corresponding to the Austwell, Mosquito 
Point, Port O'Connor and Pass Cavallo S. W. quadrangles (Figures 1 
and 2). Elevations of 5 feet or less above Mean Sea Level are generally 
encountered, with the exception of the Matagorda Island dunes and 
barrier flat, and the rangeland northwest of Port O'Connor. 

The central feature of the Austwell quadrangle, the Guadalupe 
River delta consists of deltaic muds and sands and associated 
levee and crevasse splay deposits of mud, silt and sand (McGowen, 

J. M., ec al, 1976). Marshes and grasslands with scattered shrubs 
are present. The river and its distributaries are bordered by 
thin belts of woodland which are detectable on LANDSAT imagery 
by the human interpreter, more on the basis of form than radiance. 

The areas at the junction of San Antonio and Espiritu Santo 
Bays, shown on the Mosquito Point quadrangle, include salt marshes, 
tidal flats with sparse vegetation and dredge spoil from the 
Victoria Barge Canal and Gulf Intracoastal Waterway. Vegetated 
spoil is distinguishable on the basis of form and position, in 
addition to radiance, as is also the case for beaches vs. barren spoil . 


6 



o g 

e| 






/ \ 

/ -iv ' ■ 


■ X' 


x'.-V- ’'A,,-' -■/ X' 


TEST SITE 3 






v'^X 


/ 

A 






/ 


/ 


X-. 


^ >X' 

J ) ■ ■ ■ V' !. 


.A"' 


"V *‘ - r.. 


nxM-" CrX.A XX)4 o.,:4x--. ^ i '-■' (■:■ 

-r— It r V P 'mt” ' . ,::• . 

'4 


- T. 


f > .•' ■' 


I 'vrK j 

A. 'lX 

"VJ '''>>., 

r^ AsiA' A. 


L 


■■-• ) ‘ 

/ i, 




I I 

I i 


i I A'AX:::::-- 
I i ■■- ' 


\\v "i T , X' 

'U A- XxA-:X;'^Xr^X 
y},“ T” V /■' A' ■ 

'X; ' V-A ■ X'.-XX v'- / 


■\ L.V ■ . 


II 

.'•t\ 

/ i; 


'lX-tA-I m-'"' 7^ 

/ r| \ VX V \ \ A- - ^ 

A V\VX^ A, ,:, i.\ 


'"AAA^AA' ' ■ 

' vX AfX- V^'- X ■■ 


\ 


;i:i - ' •■ ■> i ■ 


' : \ V r/ ■ \ .' 

\. ■- \ \ i, ■ K . X '■ 


l .-x 

', \ ' 

, (■■-'■ 
i ■’i 

\ I'- 

IvA>a 

. \ 

\ 

V _J. 

if 

A" a t' 

’•' 1 

r' ' 


Ij 

f 


L -V'' 


29 MARCH 1974 
SCENE 1614-16261 


/ 
r; 

'■■, XCA A 
■ ■•■■■X. I 


Figure 2. 



STATUTE MILES 

4 8 12 

-J 1 I 

SCALE 


Line boundary map showing each of 
the distinctive units for which a 
classification decision was made in 
site 3, 



The flood-tidal delta marshes of the Pass Cavallo area are the 


most extensive salt marshes in test site 3 and are the most simi- 
lar in physiography to the well known low marshes of the southeastern 
United States, The dominant species are Spartina alterni flora , 

Batis maritima and Salicornia spp, in the low areas most subject 
to tidal inundation. Black mangrove ( Avicennia germinans ) is also 
present but not in stands extensive enough to be detected on 
LANDSAT imagery. 

Cropland in site 3 is concentrated over Pleistocene inter- 
distributary mud with distributary silts and sands (McGowen, J. H. 
et al., 1976). Pleistocene marine deltaic sand, barrier strand- 
plain sands and sheet sands are not favorable for crops and 
therefore such areas are used as range-pasture lands. The latter 
include Ingleside sands southwest of Port O'Connor and south- 
southwest of Austwell, which appear mottled on LANDSAT imagery due 
to accumulation of mud and organic material in isolated low de- 
press ions. 

3.1.2 ADP Reclassification of Test Site 3 

During this quarter, test site 3 was reclassified for two 
reasons: (1) the capability of producing registered printer 

maps Is now available, and (2) a uniform, automatic data processing 
(ADP) classification for all areas within the site is desirable 
for correlation purposes. Previously, each area had been class- 
ified separately. The new classification used ISOCLS to produce 
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the class statistics and ELLTAB to classify the data. 

The use of ISOCLS to develop the class statistics is the most 
critical part of the classification procedure, since the statis- 
tics determine the sensitivity of the classification. The initial 
ISOCLS run was made using the following parameters: 

(1) one sixth of the area, every third line and every second 
sample, was clustered; 

(2) the minimum number of pixels per class, NMIN, was set 
at of the pixels clustered; 

(3) the number of iterations, ISTOP, was set at 10; 

(4) the maximum number of clusters, MAXCLS, was set at 30; 

(5) the maximum allowable standard deviation, STDMAX, was 
set at 3.0; and 

(6) clusters separated by a distance of less than DLMIN=2.0 
were combined. 

With these input parameters, ISOCLS produced 19 classes, three of 
which were small, scattered classes with large standard deviations. 
Problems working with such classes are discussed in Section 2.0. 

The statistics for the classes were used by ELLTAB to classify 
the areas of interest within the site. A recently modified version 
of DAM, the Detection and Mapping Package (for water) developed 



by NASA, Johnson Space Center (JSC), was then used to produce 
printer maps scaled and registered to 7 1/2 minute U.S.G.S. quad- 
rangle sheets (Appendix A), The changes to DAM consisted of the 
modification of one program, MAPRNT, and the replacement of one 
subroutine, READ2N, by another, RDCLAS. Some of the 19 classes 
were combined in printing the registered maps in order to facili- 
tate correlation with categories on other map data, such as the BEG 
Coastal Atlas series. 

3.1.3 Correlation of ADP Spectral Classes with Image Interpret 

tation Results 

Enlargement of portions of the 1:125,000 image interpre- 
tation line boundary map (Figure 2) to the 1:24,000 scale of the 
automatic data processing (ADP) display has greatly aided correla- 
tion of the two products. The enlargement is made from the line 
boundary map scribe sheet, that is projected on a photographic matte 
film which will accept ink or pencil annotation. The alphabetic 
abbreviation for each class on the classification scheme is added 
for each of the delineated areas, and the sheet is then overlain on 
the quadrangle-registered, ADP display for correlation purposes 
(Figure 3). This correlation techniques was conceived by Mr. Bill 
Hupp (TWDB), and the overlay produced on equipment of the Texas 
State Highway Department. 

Registration between the film overlay, the Austwell quad- 
rangle, and the registered ADP display was good, though not perfect. 
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Figure 3. Correlation on part of Austwell quadrangle using an enlarge- 
ment of the image interpretation line boundary map overlayed 
on the scaled and registered ADP display. The large numbers 
refer to results of the correlation process mentioned in 
Section 3.1.3. Hand written symbols on this working copy 
refer to the Glassification scheme: eg., Wlm=Low marsh, 

Whm-high marsh, A=agri cultural land, G=grassland, W0m=mixed 
woodland, WAmt=moderately turbid water, WAt=highly turbid water, 
Bu=undifferentiated barren. 


n 



considering the production of the film overlay from unrectified 
LANDSAT data that was enlarged* first from 1:1*000,000 to 1:125,000, 
and then from 1:125,000 to 1:24,000. The large numbers on Figure 3 
refer to the following observations that have resulted from the 
correlation process: 

1) A narrow strip of fluvial woodlands, verified from 
aerial photography (Mission 300) and field work was 
mapped along the Guadalupe River on the image inter- 
pretation, but not on the ADP product. This class 
was omitted from the ADP results, possibly because 
the radiance was not distinctive enough and the sample 
of woodland too small, relative to an adjoining grass- 
land area. 

2) A very high reflectance, barren area in a Grass- 
land area was classified as undifferentiated 
barren land on the image interpretation and was 
verified from aerial photography, but was not 
delineated as a separate class on one of the earlier 
versions of the ADP product. This barren area was 
resolved with further ADP work. 

3) Two ADP classes generally corresponded in 

location (with some other ADP classes scattered 
throughout) to the image interpretation classes 
of topographically low, wet marsh and topographically 
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higher, less wet marsh which grades into range- 
pasture grasslands. This general correspondence of 
one ADP class per image Interpretation marsh class 
may make possible the discrimination of the two 
marsh types if this distinction is verified by field 
examination. 

4) Correspondence between image interpretation and ADP 
analysis of fallow agricultural fields was ex- 
cellent, since the bare mud substrate on the 
29 March 1974 scene provided high radiance contrast 
with surrounding range-pasture grasslands. 

Further qualitative comparisons of this type will be used to build 
a table of the types of differences and similarities between 
the ADP and image interpretation products revealed by detailed 
correlation. 

3,1.4 Field Verification Approaches 

Biological field verification is conducted as required by 
the image interpreter of the project in order to assist in cor- 
relating ADP products with imagery and also to document the marsh 
vegetation represented in the ADP and image interpretation marsh 
classes. Vegetation sampling sites are chosen within spectrally 
uniform areas on the ADP and image interpretation results, and 
then located on U.S.G.S. 7 1/2 minute topographic maps. The 
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field approach currently being performed uses the point Intercept 
method (Appendix B). 

Most sampling sites selected during this reporting period 
are normally accessible only by small boat. Due to unusually 
low tides during February, these areas have been inaccessible 
even to the smallest TPWD outboard skiffs. Thus, at this time, 
only a few field transects for vegetation types that were 
accessible by land have been accomplished. 

The geologic verification procedure includes extensive use 
of available mapping, aircraft photography and published reports 
as the basis for ground truth. In addition, time histories of weather 
conditions have been compiled for the week ending with the date of 
the LANDSAT image, as an aid in interpreting each scene. Figure 4 
shows wind data for the 22-30 March 74 period which includes the 
29 March 74 scene of test site 3 and is based on U. S. Weather 
Service data taken at Victoria, Texas, As this station is some 35 
miles inland from the test site the precise velocities will not 
apply to the coastal area* but the general onshore or offshore 
direction and relative duration are probably similar. Precipita- 
tion data is also being compiled and tide gauge data has been 
obtained from the li. S, Army Corps of Engineers, Galveston District. 

The large scale photography of Missions 300 and 325, as well 
as the smaller scale photography for the entire Texas coastal zone, 
have been very useful in interpreting questionable features on the 
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Figure 4. wind data prior to the 29 March 1974 image shows 

that wind strengths were low, and therefore, water 
levels would not be wind influenced. Effective winds 
are those above 10.4 knots (12 mph), and S=(V-10.4 
knots)2(d) where V=observed velocity and d=wind dur- 
ation (Price, 1975). Data are plotted at 3-hour 
intervals, and shaded areas indicate periods of 
effective winds and the associated peak strength. 
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LANDSAT imagery. When this photography is utilized with the 
Environments and Biologic Assemblages and Environmental Geology 
sheets of the Coastal Atlas (McGowen, J. H.» et al., 1976), an 
excellent natural resource base of information exists for test site 
3. It is anticipated that one additional field trip will be made 
for geologic verification following the interpretation of all 
scenes for this test site. 

3.2 Examination of Test Site 2 (West Galveston Bay Area) 

3.2.1 Preliminary Site Description and Image Interpretation 
Results 

After examination of topographic maps. Environmental 
Geologic Coastal Atlas maps (McGowen, J. H., et al., 1976) and 
NASA aircraft photography covering the Freeport-West Galveston 
Bay area (Site 2), a visit to test site 2 was made on 27-28 
January 1976. The objective was to gain some familiarity with 

patterns of land use, geologic materials and the distribution 
and general types of vegetation prior to LANDSAT image inter- 
pretation. Test site 2, located about 60 miles northeast from test 
site 3, is probably the most diverse among all the test sites in 
that it includes: 

(T) the urban areas of Freeport, Lake Jackson and Clute, as 

well as numerous smaller developed areas; 

(2) major chemical and petroleum-production related indus- 
trial sites along with numerous smaller sites for 
industries serving the petrochemical complexes; 
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(3) tha Brazosport shipping facilities; 

(4) many dredged channels associated vnth the area's 
commercial activities; 

(5) a system of hurricane protection dikes resulting in 
non-natural contrasts in vegetation type over short 
distances; 

(6) andj extensive v/etlands associated with a strand- 
plain type of shoreline developed on a Holocene deltaic 
headland (McGowan and Scott, 1975). 

This complexity has resulted in the large number of individual 
areas delineated on the line boundary map (Figure 5), interpreted 
from scene 1289 - 16261, 8 May 1973. The presence of many active 
and abandoned stream channels, ponds, oxbow lakes and low marsh 
areas resulted in the delineation of numerous small units on the 
basis of form and texture. 

3.2.2 Preliminary ADP Classification 

ADP classification of test site 2 (Scene 1289-16261, 8 May 
1973) was begun during this quarter, and followed the procedures 
described for test site 3 in the preceding progress report. DAM 
was used to classify water and to register the classified results 
to U.S.G.S. 7.5 minute topographic maps. ISOCLS was used to pro- 
duce a uniform set of statistics for the site. Again, the problem 
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Figure 5. A line boundary map of site 2 showing 

the complex active and abandoned drainage 
patterns associated with this strandplain 
shoreline. 



previously noted, that of small -member, scattered classes with 
large standard deviations, was encountered. ELLTAB was used 
to generate a preliminary classification of the area. Land 
features and water features were classified in separate runs 
when the combined look-up table exceeded the storage alio ted 
to it. Registered classification maps are currently being 
produced for test site 2 that contain some areas of unclassified 
data. As soon as these difficulties with the poor statistics 
have been resolved, another classification at the site will be 
produced for correlation with the image interpretation results. 

3,3 Interim Status of Monitoring System Design 

3.3.1 Data Acquisition 

The current status of LANDSAT data acquisition to support 
this investigation is reflected in Appendix C. The primary em- 
phasis was placed on acquiring CCT's for sites 4 and 5 based on a 
review of the corresponding imagery which was previously ordered. 
In addition, all required imagery and tapes were ordered for 
LANDSAT scene 2375 - 16112, (1 February, 1976) covering site 2. 
This will provide a good winter scene of the site using the most 
current coverage which can be received within the time constraints 
of the project. 

The "Suggested Purchase Package for LANDSAT Imagery Scenes" 
(Table 1) provides the necessary imagery to satisfy the imagery 
interpretation aspects of this investigation. In addition to 
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Table 1 


Suggested Purchase Package For 
LANDSAT Imagery Scenes 



Item 

Scale 

Band 

Cost 

1) 

color transparency 

1:1,000.000 

composite 

$12. 00 

2) 

B & W print 

1:250,000 

5 

15.00 




7 

15. 00 

3) 

B & "W positive 

1:1,000,000 

4 

5.00 


transparency 


5 

6. 00 




6 

5. 00 




7 

5. 00 

4) 

B & W negative 

1:1,000,000 

5 

6. 00 


transparency 


7 

6. 00 


Cost if color master available ...... 


. .$74.00 


Additional cost of color master ...... 


. $50.00 


Possible cost if color master needed $124.00 
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these requirements, the ADP classification effort requires the 
Computer Compatible Tapes (CCT) and a black-and-white positive 
transparency of Band 7 to support the Detection and Mapping (DAM) 
classification of land vs. water. 

Experience gained from previous orders of LANDSAT materials 
placed with the EROS Data Center (EDC) indicate the following de- 
livery schedules; Black and white imagery of the 4 bands and the 
tCT‘s require 4-6 weeks from the date of acquisition by the initial 

* N 

processing site (Goddard) for delivery of a ''master" copy to EDC. 
Once these copies are available at EDC, it requires 2-3 weeks for 
production of black and white or color composite imagery products 
if the color composite "master" has already been produced. Pro- 
duction of the color "master" will require an additional 2-3 days. 

Delivery Of the CCT's requires 7-10 days. Shipping time will 
probably add 2-3 days to each delivery. These time frames quoted 
above are generally maximums and in most cases delivery will be 
made in less time. 

3,3.2 Interim Status of ADP and Image Interpretation Techniques 
A revision of the ADP classification schedule presented 
in the last quarterly report is included in Table 2. Except for 
some changes in the parameters used In step 9 for developing 
class statistics using ISOCLS, only minor revisions in the ADP 
Classification schedule are anticipated during the next quarter. 
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Table 2 


LANDSAT ADP Classification Schedule 

Scene ID; 

Test Site: 

Description: 

1. Select LANDSAT scene and determine data tapes ID number. 

2. Examine available imagery. 

3. Estimate scan line and sample numbers for the areas of interest. 

4. Merge data tapes if necessary. 

5. Generate grayscale maps of the area. (GRAYMAP) 

6. Obtain meteorological data (precip, wind, tides). 

7. Locate control points for the scene for registration and scaling, 

8. Classify water using DAM. 

9. Cluster all training areas within the scene (ISOCLS). 

TO. Examine class statistics. 

11. Refine a training class if indicated by step 10. 

12. Use class statistics to build the look-up table (ELLTAB). 

13. Classify the area (ELLTAB). 

14. Register and display the classified results (REGISTER). 

15. Outline or color code homogeneous areas. 

16. Examine the classification map. 

17. Stop if satisfied with the results. 

18. Retrain on unclassified or poorly separated areas (ISOCLS). 

19. Go to step 12. 
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The techniques thus far developed for LANDSAT Image Inter- 
pretation at a 1:125,000 scale have formed the basis for the schedule 
of analysis of each scene outlined in Table 3. Note that a general 
acquaintance (steps 1 and 2) with the area is considered necessary 
before the line boundary delineation and area classification 
procedures are undertaken. Data review during step 1 is not as 
detailed as in step 5 during the verification procedure. 

The accuracy evaluation procedure (step 8) will involve a 
qualitative analysis of the differences between LANDSAT image 
interpretations, ADP products and a third source of information. 

This latter will primarily be a combination of the NASA large scale 
aircraft photography of each test site and maps from the Environ- 
mental Geologic Coastal Atlas series (University of Texas, Bureau 
of Economic Geology, 1972-). Results will include a table of the 
types of differences between LANDSAT and larger scale data and the 
frequency with which they were encountered. The table will include 
each of the categories in the land use classification scheme for 
i mage i n terp re t a ti on . 

3.3.3 Interim Display Products 

At this time the products that we are considering to 
display the image interpretation and ADP classification efforts 
will be of two types: 

1) Printer maps scaled and registered to US6S 71/2 minute 
topographic maps (1:24,000 scale) showing the ADP 
classification. These maps must be hand colored. 
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Table 3. 

'Schedule For Image Interpretation Analysis 


1. Review aerial photography. Coastal Atlas Maps, and published 
tide and weather data for test site and image date. 

2. Take a preliminary field trip to become generally acquainted with 
test site (may include oblique aerial photography). 

3. Complete line boundary map of test site area. 

4. Classify features delineated according to the modified Anderson 
system. 

5. Study supportive data in detail, review results, field check and 
correlate with biological verification. 

6. Document results for the scene, especially problems and unique 
aspects of the imagery. 

7. Produce corrected image interpretation at 1:125,000 scale and 
overlays of selected quadrangles at 1:24,000 scale. 

8. Qualitative analysis of the classification products to evaluate 
accuracy, utilizing aerial photography, Coastal Atlas maps and 
field results. 

9. Evaluate format and content of resulting map. 

10. Evaluate image interpretation of this scene in conjunction with other 
scenes for the same test site. 
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taped together, and trimmed to be used. In addi- 
tion, the information on these maps must be further 
interpreted to make them completely compatible with 
the classification scheme that we have developed. 

2) The image interpretation display product (1:125,000 
scale) will be reproduced on a stable white plastic 
material, and hand colored to show the classification 
of all areas delineated on the line boundary map. 

For this product the interpretation will be based on 
the classification scheme. 

Since this investigation does not contemplate generating a 
map series, these products or interim versions similar to these 
will be all that we will attempt. 

The TNRIS however, is investigating alternative ways of 
generating hard copy color final products of ADP results, so that this 
service would be available to users of TNRIS. For example, TNRIS 
has requested that Sei scorn Delta, of Houston, Texas, generate a color 
version of the Austwell classification developed in this investigation, 
and also submit a price list for different scales and volumes of 
products following this theme. 

The products discussed above will be those compared to other 
information products of similar nature in the cost-savings analysis. 
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3.4 Program for the Next Reporting Interval 

The following activities are planned for the next 
quarter: 

1) Techniques developed for analyzing LANDSAT data 
will be applied to test site and also at least 
one additional scene in test sites 2, 3, and 5, 
will be analyzed to determine the effect of sea- 
sonal and other changes in these scenes, 

2) The optimum technique or mix of techniques for analysis 
of the LANDSAT data will be evaluated for use on the 
Texas coast. 

3) Design of the LANDSAT-based monitoring system will be 
completed for testing during the summer. 

4.0 SIGNIFICANT RESULTS 

The most significant ADP result during this quarter was the modi- 
fication of the DAM package to produce classified printouts, scaled 
and registered to U.S.G.S., 7 1/2 minute topographic maps from LARSYS- 
type classification files. With this modification, all the powerful 
scaling and registration capabilities of DAM become available for 
multi-class classification files such as those produced by LARSYS and 
ELLTAB (Section 3.1.2), 

The most significant results with respect to image interpretation 
have been the application of the mapping techniques to a new, more 
complex area (test site 2), and the refinement of an image interpretation 
procedure which should yield the best results (Section 3.3.2). 
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5.0 PUBLICATIONS 

The ADR classification of the Austwell quadrangle in site 3* generated 
by the procedure reported In this progress reports will be processed 
by Sei scorn Deltas Inc. of Houston, Texas, during this next quarter. 

Seiscom Delta, Inc. will produce a sample color photographic product 
at a reduced scale and prepare a price list so that this option will 
be available to Texas agencies who plan to access the ADP classification 
procedure for LANDSAT data developed by this investigation through 
the Texas Natural Resources Information System (TNRIS). 

6.0 RECOMMENDATIONS 
None . 

7.0 FUNDS EXPENDED 

GENERAL LAND OFFICE {6L0) 

- Labor 
Overhead 
Travel 

TOTAL EXPENDITURES REIMBURSED DURING THE 
3rd QUARTER 

BUREAU OF ECONOMIC GEOLOGY (BEG) 

Labor 

Materials & Supplies 
Equipment 
Travel 

TOTAL EXPENDITURES REIMBURSED DURING THE 
3rd QUARTER 


$ 6,909.00 
6.90 
163,04 

$ 7,078.94 

$ 0 

0 
0 
0 

$ 0 
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TEXAS PARKS & WIIDLIFE DEPARTMENT (TPMD) 


Labor 

$ 0 

Materials & Supplies 

n 

w 

Travel 

0 

TOTAL EXPENDITURES REIMBURSED DURING THE 
3rd QUARTER 

$ 0 

TEXAS WATER DEVELOPMENT BOARD (TWDB) 


Labor 

$ 2,822.00 

Computer 

400.00 

TOTAL EXPENDITURES REIMBURSED DURING THE 
3rd QUARTER 

$ 3,222.00 

CONSULTING SERVICES 


Dr. John A. Schell 

$ 678.96 


TOTAL EXPENDITURES REIMBURSED DURING THE 

3rd QUARTER $ 678.96 


CUMMULATIVE TOTAL EXPENDITURES REIM- 
BURSED DURING THE 3rd QUARTER $10,979.90 
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8.0 DATA USE AS OF FEBRUARY 29, 1976 



IMAGERY 

Account #623790 
Amount 

CCT 

Account #G B3790 
Amount 

AIRCRAFT 
Account #G W3790 
Amount 

Value of Data 
A1 1 owed 

$2,900.00 

$5,400.00 

$9,588.00 

Value Ordered 

$1,465.00 

$2,000.00 

$9,564.00 

Value Received 

$1,274.00 

$ 600.00 

$9,564.00 

BALANCE 

$1,435,00 

$3,400.00 

$ 24.00 


9.0 AIRCRAFT DATA 

The NASA aircraft photography has assumed greater importance with 
the development of the image interpretation sequence (Section 3.3.2, 

Table 3). It serves to familiarize the interpreter with the study 
area so that reasonable line boundaries can be drawn and classification 
decisions made. In an operational system the 1:120,000 photography of 
the entire Texas coastal zone (Mission 300, February, 1975) would fulfill 
this function. The NASA aerial photography also is being used in the 
verification procedure and will be the primary tool for analysis of 
classification and boundary accuracy as study of each scene is completed. 
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APPENDIX A 

SUBROUTINES MODIFIED IN DAM TO PRODUCE PRINTER MAPS 
SCALED AND REGISTERED TO U.S.G.S. 7 1/2 MINUTE TOPOGRAPHIC MAPS. 



0B0200-q2*SCAUE.ROCLAS 

1 SUnROllTiNE RDCLAsnLDISKflSDWLOflsnWItlftJXWD) 

2 C4^**TH1S subroutine reads a uarsys cuassification file 

3 C+**AND RETIIRF4S CLASSIFICATION RESULTS FOR LINE ILDtSK 

4 DIMENSION NXWn(4C300) tlDATA(lOOO) rNPTS(4) rLSTART(4) fLEN0(4) t 

5 +15TAf?T{4) ♦ IENn<4) 

6 DATA JUMP,T1SAM/0»010/ 

7 DEFIHF XCCT(5nTSK)={SDISK+NSAH-l)/NSAM 

8 IF (JlJMP.EQ.l) GO TO 10 

9 J1IMP=1 


10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 


NCCTLO=NAXtUXCCTaSDWLO~l) >1) 
NCCTHI=MlMOtXCCT(ISDWHl+l) t4) 
C**+REAO Hr.ADER RECOf^US FROM LARSYS FILES 
DO 5 MCCT=MCCTl-0,NCCTHl 
NLI!!IT=24+NCCT 
RtADFrJUNlTl MXWO(i) 

READFMIJNIT) NXWOai 
READtNllNIT) NXWD(U 
READ(NIINITV NXWO(l) 

READFNUNIT) ( NXWD { I W I=l» 10 3 
C*+*SAVE FIELD INFORMATION 
NPTSIMCCT)=NXWR<1) 
LSTART(NeCT)=NX.;Dt6) 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


LEN0(NCCT)=NX'-;r.t7) 

ISTART(NCCT)=NXi‘JD(9) 

5 IE(JO(NCCT)=NXWn(10) 

C'^**FlLL LINE WITH «NO DATA* FLAGS 
10 DO 20 I=lf4000 • 

20 NX’/D{ I) = *000000' 

C+*=KLOCATE REQUESTED DATA ON CLASSIFICATION FILES 
DO 50 NCCT=NCCTLQ»NCCTHI 
NUNIT=24+NCCT 

IF (ILD1SK,LT.L5TARHNCCT3 .0R*ILDISK.GT,LEN0CNCCT)) GO TO 50 
C**+READ A CLASSIFIED LINE 
J=NPTStWCCT) 

30 READtNUNITf£NO=50) ILlNEf CIDATACl)il=l>J3 

IF- ( JLINE.LT.ILDISK) GO TO 30 
C+*+InSERT CLASSIFIED DATA INTO NXWD 
J=0 

KK=ISTART { NCCT ) +NSAM* ( NCCT-1 ) 

LL=I END ( NCCT 3 f NSAM* ( NCCT-1 3 

DO 40 I=KKfLL 

d=J-t-l 

40 NXWD(I)=IDATA(J3 

50 CONTINUE 

RETURN 

C**>!=RErtrNO FILES AND RESET JUMP FLAG 
ENTRY RE5ETJ 
JUMP=0 

00 60 NCCT=NCCTLO»NCCTHI 
NUMIT=24+NCCT 
60 REWIND NUNIT 

RETURN 

END ■ ■ 
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DB020 0-02+SC ALE , MAPHNT 


1 

SURROUTItJE MAPRNTtKERR) 


2 

C (EriS) 


3 

C THIS SUBROUTINE REGISTERS ERTS MSS DATA FOR 


4 

C DAM.PRTCLASS. 


5 

C 


6 

c 


7 

C THIS SUBROUTINE CALLS THE FOLLOWING EXTERNAL SUBROUTINES/FUNCTIONS: 

B 

C NITHDG 


9 

C SYMTAB 


10 

C TICGFM 


11 

C RDCLAS 


12 

C 


13 

INCLUDE KOMXOT/LIST 


14 

INCLUDE KOMIDrLIST 


15 

INCLUDE KONCTL»LIST 


16 

INCLUDE KOMCCTfLlST 


17 

INCLUDE KOMMAPfLlST 


13 

include komaltvlist 


19 

INCLUDE KOMSYMrLIST 


20 

DATA IOUT/8/ D OUTPUT TAPE OF CLASSIFIED 

DATA 

21 

OltCNSIOU HXWD(40001 


22 

DIMENSION rSpSYMC112G) 


23 

OIMEfSSiCiN LlNFMT(4> 


24 

DATA LINFMT/' (lK,Il»lXrJ4rNNNAl»J4) */ 


25 

DATA JDOLAR/'S*/ 

FOR TAPE correction 

26 

INCLUDE TAPCORrLIST 0 DEFINE PROCEDURES 

27 

include N1 TAB, LIST Q DEFINE PROCEDURE TO 

COMPUTE ALT PRINT UNIT 

2S 

INCLUDE DIGITStLIST (3 DEFINE PROCEDURES 

FOR DIGIT EXTRACTION 

29 

C 


30 

c 


31 

C generate TABULAR DATA 


32 

C 


33 

NITLOm 


34 

NITHI=0 


35 

INCLUDE NITROTfLIST 


36 

NIT=:C1 


37 

NUNITrNTADtNITJ 


38 

CALL NITHDG(MIT,NUNITJ 


39 

CALL SYMTABtNIT,NUNlT) 


40 

CALL TIC6EN(KERR»NIT,NUNIT) 


41 

C 


42 

C 


43 

C INITIALIZE WINDOW 


44 

C 


45 

rLPWIN=WLPMlN 


46 

rLPWAX=,-.'LPMAX 


47 

ISPWIN=WSPMIN 


43 

ISPWAX=isSPMAX 


49 

ILPTIC=C) 


50 

ISPTIC=0 


51 

LVLTIC=» 


52 

c 


53 

C 


54 

C BREAK WINDOW INTO SUB-WINDOWS 


55 

c ■■ ■ ■ ■ 


56 

ISPMOO=MOO( (ISPWAX-I5PWIN) ,1241+1 
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HNOFILE I OUT 

172 



CALL RESETJ 

173 



RETURN 

174 

C 



175 

C 



176 

C 



177 

C 



178 

C 



179 

C 



160 

C 



181 



SUBROUTINE GETIC 

102 



MTICKSNTICK+I 

163 



ILPTIC=FLD{OOrlO»LOCTIC(NTICK) ) 

104 



ISPTIC=FLD ( 18 f 12 1 LOCTIC (NTICK ) 5 

165 



FLD{0»6fLVLTIC)=FLOt30»Q6iL0CTIC(NTICK) ) 

166 



RETURN 

167 

c 



188 

c 



189 

c 



190 

c 



191 

c 



192 

c 



193 



SLtBROUTINE SA^BCL 

194 


99 

FORMAT (7X* 12411) 

195 



ISPNLOnSPWLO 

196 



DO 990 NIT2NITL0,MITHI 

197 



NUM1T=UTA9(NIT) 

198 



ISPNHI=MIi;0( aSPNLO+123) tISPWHl) 

199 



DO 1000 i=ispmlo,ispnhi 

200 

1000 

IF (I-ISPWIN+l.GT.O) 1SP5YM(I-I5PWIN+I)=l/H)0n 

201 



VJPITE (MUIJITj 99) ( ISPSYM 1 1-ISPWTN+l ) t I-ISPNLOf ISPMHI ) 

202 



DO 100 r=lSPNLO»ISPNHl 

203 


100 

IF tr-ISPWIN+l.GT.O) ISPSYM(I-ISPWlN+l} = n-1000+(l/10001 )/100 

204 



WRITE {MUNI T» 99) t ISPSYM { I- ISPWIM+1) f I=ISPNLO» ISPNHI ) 

205 



DO 10 I=lSPNLOrISPNHI 

206 


10 

IF (I-ISP'VIN+I.GT.O) I5P5YM(I‘-ISPWIN+1) = (I-100*(I/100))/10 

207 



WRITE{riUNITf99) (ISPSYMU-I5PWIN+1) »I=ISPNLO» ISPNHI) 

206 



DO 1 I^ISPNLOrlSPNHI 

209 


1 

IF { I-ISPWIN+l.GT.O) ISPSYM(I-1SPWIN+1)=I-10*(I/10) 

210 



WRITE(NUNIT»99) (ISPSYM£I-ISPWIN+1) t I=ISPNUO* ISPNHI) 

211 


998 

ISPNL0=ISPMHI+1 

212 



DO 999 I=ISPWLO.ISPWHI 

213 


999 

IF (I-ISPWIN+1.GT,0) I5P5YM{I-ISPWIN+1)=0 

214 



RETURN 

715 

c 



216 

c 



217 

c 



218 

c 



219 

c 



220 

c 



221 



SUBROUTINE LINOUT 

222 

c 



223 

c 

LOOK UP COUNT SYMBOLS 

224 

c 



225 



WRITEtlOUT) { ISPSYM ( I-XSPWlN+1 ) »I=XSPWLO,ISPV'JHl) 

226 



DO 150 I=ISPWLOrISPWHI 

227 



ISPCNT=MINO (ISPSYM (I-ISPWIN+l ) t (KSYMSZ-1) ) 
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116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 
136 
139 


C 

c 

C SCAN DENSITY PIXELS 
C 

320 ISDISK=ISDISK+1 
fiWD=NWn+l 

325 IF(ISDTSK.GT.ISDWHI) GO TO 400 
GO TO 310 
C 
C 

C REGISTER FIRST *NO DATA' PIXEL 
C 

350 AGCOt^rXASCOR (FLOAT! ISniSK) ) 

I SP 1=CP { 4 ) +ALCOR+CP { 5 ) *ASC0R+CP ( 6 ) 
IF(lSniSK.GT.O) GO TO 360 
iSniSKrl 
NWD=1 
C 

c 

C SCAN »NO DATA* PIXELS 
C 

360 lSDISK=ISniSK+l 
NWD=NWt3+l 

365 IFCISOISK.GT.ISDWHI) GO TO 380 

IF (NXWD(NWD) .ME. *000000*) GO TO 380 
GO TO 360 


140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 


C 

C 

C REGISTER STRING OF *NO DATA* PIXELS 
C 

330 ASCORrXASCOR {FLOAT ( ISOISK-1 ) ) 13 LAST *NO DATA* PIXEL 

1SP2=CP ( 4 ) *ALCOR+CP ( b) *ASCOH+CP 1 6 ) 

DO 385 1SPRNT=ISP1»ISP2 

385 IF (ISPRNT-ISPWIN+I.GT.O) ISPSYMtISPRNT-ISPWIN+l)=+999999 
IFdSDlSR.GT.ISOwHll GO TO 400 
GO TO 310 
C 

G . 

C INCREMENT DISK LINE AND WRITE PRINT LINE 
C 

400 ILDISK=1LD1SK+1 

ALCORsILDlSK IDFUTURE CORRECTION 
NLPRNT=CP ( 1 ) *ALC0R+CP ( 3 ) 

IF(NLPRNT.GT.ILPRNT) CALL LINOUT 
IF(NLPRNT.GT.ILPwAX) GO TO 500 
GO TO 200 
C 
C 

C FOOT SUB-WINDOWS 
C 

500 CALL SAM5CL 
C 

c 

aOO CONTINUE 

WRITEUJUNIT.aOS: 

805 F0RMATt*0*/6X.*#»SEE UNIT 0 FOR legend**') 

C 


J-.IGtNAL PAOR IS Prv-:-. 





5a 

59 

oO 

61 

62 

63 

64 

65 

66 
67 
6S 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
BO 
ai 
82 

83 

84 

85 

86 
67 
08 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
113 


FLn(r;6,6»LIUF'-lT(3) , rj- I.D ( 0 0 r 6 r JMtJijS C ISPMOf') ) 

FLD < 1 a » 6 » L I WF VT 1 3 ) ) =FLn ( 0 0 , o ► JTCNS C 1 5PMOD ) ) 

FLO{15»6*LINF‘‘iT(i) )=FLDt00»6* JONLSCISPMOO) 1 

D9 800 NlTLO=l»niTMAXfMALTM 

|jlTHI=-'Itm C (NITLO+N'ALTM-1 J ,MITfJ'AX) 

ir'CLPDF SflTKOT^LIST 

ALC OH =PC (1 ) * I LPW I N+PC ( 3 ) 

ILOISK=ALCOK 0 FUTURE CORRECTION 

alcorsildisk q future correction 

I LPRNT=CP < 1 ) +AUCOR+CP ( 3 ) 

ISPWLO=ISPWIN+ 124=*' CNITLO-1 ) 

ISPWHI=MIN0 t ( T5PWIN+124WITHI-1 ) r ISPWAX) 

NTICK=0 
CALL SETIC 
C 

c 

C HEAD SUB-WINDOWS 
C 

00 140 NIT=NITLOrNITHI 

NUNIT=NTAB(NIT5 

CALL NlTUDGtNlT.NUNlT) 

140 CONTINUE 
CALL SAMSCL 
C 

c 

c COMPUTE FIRST/LAST DENSITY SAMPLES 
C 

200 ASCORrPC { 4 ) *ILPPNT+PC (5) *ISPWLO+PC (6) 

ISDWLO=X ASnSK ( ASCOR ) 

ASCOR=PC ( 4 ) * ILPRNT+PC 1 5 ) * ( ISPWHI+1 J +PC (6) 

1SDV)HI=XASDSK ( ASCOR ) +1 . 0 
C 
C 

C READ DENSITY LINE 
C 

CALU ROCLAS ( ILDISK » ISOWLOr ISDWHI rNXWD) 

C 

c 

G LOCATE FIRST DENSITY PIXEL 
C 

ISDISKrlSDWuO 
NwDLO=tSDWLO 
NWO=f!WDLO 
NWDHI= ISDWHI 

IFdSDISK.LT.l) GO TO 350 
C 
C 

C SCREEN PIXEL DENSITY 
C 

310 IF (NXWDCNWDl.EQ.’OOOOOOt) GO TO 350 
C 
c 

C REGISTER/COUNT SCREENED PJXELS 
C 

ASC0R=XASC0H (FLOAT (ISO rSK)) 

I SPRNTSSCP { 4 ) ALCOR+CP (5 > * ASCOR+CP (6 ) 

IF USPRNT-ISPWIN+l.GT.O) ISPSYMMSPKNT-ISPWIMn)aNXWD(NWD) 



223 

229 

230 

231 

232 

233 

234 

235 

236 

237 
233 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 
26/ 
268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 
283 
204 


150 IF tI-ISPWlW+l,GT.03 ISPSYHCI-ISPWlN+l>=KSYM(lSPCNT+lJ 
GO To 300 
C 
C 

C tJLANK OUT LirJE SKIPPED FOR SCALING 
C 

200 DO 250 I=ISPliLO. ISPWHT 

250 IF (I-ISP'7IN+1.GT.0J 1SPSYM(I-ISPWIN+1) = * 

C 

c 

C INSERT TICK MARKS 
C 

3QD IF(ILPTIC.GT.ILPRNT) RO TO 400 

IFtlSPSYfinSPTlC-ISPWlM+D.EQ.’ n GO TO 330 
IFtI5PSYMnSPTIC-ISPwIN+l) .NE.* ; • ) GO TO 350 
IF(I5PSYM{ISPTIC-ISPWIN).EG.» :« ) 

1 ISPSYM{ISPTIC-ISPWIN)=* * 

IF ( ISPSYMUSPTIC-ISPwIN-)-2) .eo. 'in 
1 15PSYJ^{I5PTIC-lSPWlN+2) = * • 

330 IF (ISPTIC-ISPWIN+I.GT.O) ISPSYMUSPTtC-lSPWlN+l)=LVLTIC 
350 CALL GETIC 
GO TO 300 
400 CONTINUE 
C 

c 

C WRITE SUU-WINDOW PRINT LIME 
C 

ISPNLO=ISPWLO 

DO 540 NIT=NITLO,NITHI 

NUNIT=NTA.r!{NIT) 

ISPNMI=MIN.0 C { I5PNL0+125) r rSPWHl ) 
IFC(NlT,EQ.NITMAX),AND.aSPMOD.LT.122)) GO TO 530 
IvRITFJNUNITfSaO) NlT.lLPRNTr 
1 (ISP5YM(I-I5PWIN+l>»I=I5PNLOfISPNHI3 fJDOLAR 
520 FORMAT<lXtIlflXf J4f 125Ali 
GO TO 540 

530 WRITE(NUNITrLINFMT) NITtILPRNT» 

1 (I5P5YM(I-ISPWIN+l)f I=I5PNL0»ISPNH1J fILPRNT 
540 I5PNL0=1GPNHI+1 

IF(KSY[!lT,EQ.m GO TO 700 
C 

c 

C OVERPRINT SYMBOLS 
C 

DO 660 K[3rT=L)6»KSY8lTf6 
DO 610 I = ISH.iLOplSPWHI 
610 IF (I-TSP^ilN+l.&T.O) 

1 FLD C 00 , 6 , ISPSYV {I-lSPWlN+1 ) ) =FLD {KBXTF6r ISPSYMI I-ISPWlN+1 ) ) 
ISPNLO=ISPWLO 
DO 640 NIT=NITLOfNlTHI 
NUNIT=NTA!UNIT) 

1SPNMT=MIM0 t ( lSPNLO+123) f ISPV’iHD 
WRITC{NUNIT»62n) 

1 tlSPSYM{I-ISPWlN+l)fI=:ISPNLO#ISPNHI) 

620 F0RMATt»+'»6XFl24Al) 

640 ISPNLO=ISPMHI+l 
660 CONTINUE 
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205 C 

2Qb C 

287 C INCREMENT PRINT LINE 

288 C 

209 700 ILPRNT=ILPRNT+1 

290 irtNLPRNT.GT.TLPRNT) GO TO 200 

291 C 

292 C 

293 C REINITIALIZE LINE 

294 C 

295 DO 750 isISPWLOfISPWHl 

296 750 IF (I-ISPWIN+l.GT.O) ISPSYM( I-1SPWIN+1)=0 

297 RETURN 

298 END 


fiJPRTpS scale.rdclas 


OP POOR QUAIIT^I 





APPENDIX B 


VEGETATION SAMPLING FOR LANDSAT DATA 



Vegetation Sampling for LAHDSAT Data 


prepared by Larry Lodwick 
Texas Parks and Wildlife Department 


INTRODUCTION 


In an effort to obtain quantitative ground data on the plant cotnaunlties 
as defined by LANDSAT Telemetry » a sampling system which will allow a quick, 
yet quantitative analysis of the communities needs to be developed. The pro- 
cedure should be adaptable to the various communities present, from mud flats, 
with only sparse vegetation cover to salt marshes, being predominately grasses 
or grass-llke plants, to woodlands. 

The most efficient sampling system for large areas would be the point inter- 
cept measurement of cover. This has several advantages in that 1) it gives 
an indication of biomass (especially if height of the vegetation is known) ; 

2) it can be used for all growth forms, from bryophytes to tree canopies; and 

3) it can be adapted for the size of the community to be sampled (i.e., points 
farther apart for larger vegetation types) . Although widely spaced points tend 
to reduce the precision of the measurements, it does serve as more rapid measure- 
ment than other sampling techniques (Mueller-Dombois and Ellenberg, 1974) . 

From the data obtained by the cover method, it is possible to designate 
plant associations which could be related to the various images interpreted 
by LANDSAT. 


MATERIALS AND METHODS 


An advantage to sampling cover as opposed to other parameters is the small 
amount of equipment required for field sampling. The required materials con- 
sist of a tape measure with a minimum length of 25 meters (or 25 yards) , a meter 
(or yard) stick, tally sheets (figure 1) for the sampling data and a plant press 
for collecting those plants which the investigator is unfamiliar with. 

The method for data collection is as follows : 

1) Prior to going into the field the sampling site should be located 
using the LANDSAT printout to determine the approximate center of 
the vegetation type to be sampled. This point should then be lo- 
cated on a topographic map (D.S.G.S.) or low altitude photograph. 

2) Using the topographic map or photograph, locate the sampling site 
on the ground, setting a stake at the center point. The tape 
measure should then be extended first to the north, then south, 
east, and west of the center point (preferably with the use of a 
compass) to a distance of 25 meters or yards (figure 2). The pur- 
pose of determining the location and direction of the transects 
prior to beginning of sampling is to reduce the bias of the inves- 
tigator. 


So - io 



3) At 25 evenly spaced points in each of the four directions from the 
center point, preferably at 1 meter (or yard) Intervals, all species 
directly above or below the points should be identified and its 
height, measured with the meter (or yard) stick (which, in the case 
of trees, may be estimated), recorded on the tally sheet (figure 3 
and 3a). Bare soil, without vegetation, should also be recorded 
and treated as a species. This will enable one to assess the bare 
ground (mud flats, dredge spoil, etc.). Those plant species in 
which the investigator is unfamiliar should be collected, pressed, 
and sent in for identification. Preferably the collection should 
consist of two or three individuals pressed separately with flowers 
(or fruits) and roots. These should be pressed, aeee rd - lng - to -the 
i5fftgttgb iujis ~~' l n : appendls r-A:. Information as to the sampling site, 
soil type, soil wetness (tidal marsh, dry uplands, standing water, 
etc.) should be recorded. The unknown plants should be numbered 
and the number listed on the tally sheet in place of the name. 

After the plant is identified, the number should be replaced by 
the correct name, 

A) The information on the tally sheet should Include the name of the 
investigator, to the transect line number (as related to the map), 
the bearing of the line (north, south, east, or west), the date, 
and the amount of inundation at the time of the sampling (i.e., 
dry land, mud, standing water, etc,). 

5) One copy of each field sheet should then be sent to Austin for 
evaluation and analysis of the plant associations, 

6) After several sites have been investigated, any problems encoun- 
tered need to be discussed to determine what changes in the pro- 
cedures might be made to alleviate the problems. 


Reference 


Mueller-Dombois, Dieter, and Ellenberg, Heinz, 1974, Aims 
and Methods of Vegetation Ecology : John Wiley and Sons, 

New York, 547 p. 
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Figure 1. Sample tally sheet for recording ground cover. 
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Figure 2. Sample area for the measure of the point intercept 

method. After reaching a predetermined sampling site, 
select the center point and with the use of a compass, 
record those species which occur at 25 regular inter- 
vals (preferably one meter intervals) directly north, 
south, east, and west of the center point. 
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TAXOOlUM 

DISTICHUM 


Figure 3. A schematic representation of a wetland vegetation 
type containing seven species. 
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Figure 3a. A tallv sheet with the data for the 
association shown in Figure 3. 




LANDSAT COVERAC" 0?* THE TEST SITES 2, 3, 4, 5 
FOR LANOKAT INVESTIGATION 1^23790 



APPENDIX C 


LANDSAT COVERAGE OF THE TEST SITES 2, 3, 4, 5 


ACQUISITION 

STATUS 

FOR LANDSAT INVESTIGATION 
SCENE ID 

#23790 

DATE 

CLOUD 

COVER 

QUALITY 

2 

Test Site 2: 

Summer: June - August 

1037 - 16251 

08/29/72 

Z0% 

8888 


1343 - 16253 

07/01/73 

20% 

8888 


1361 - 16252 

07/19/73 

20% 

8888 

2, 4 

1703 - 16175 

06/26/74 

10% 

8858 

2 

Fall: Sept. - Nov. 
1073 - 16251 

10/04/72 

30% 

8888 

2 

Winter: Dec. - Feb. 

1217 - 16261 

02/25/73 

20% 

8888 

2 

1505 - 16230 

12/10/73 

00% 

2822 

1 

1901 - 16110 

01/10/75 

10% 

8808 

1, 3 

2375 - 16112 

02/01/76 

00% 


2 

1576-16152 

02/19/74 

00% 

8888 


Spring: Mar. - May 
1253 - T6262 

04/02/73 

20% 

8888 

2, 4 

1289 - 16261 

05/08/73 

. 00% ^ 

8888 

2 

2051 - 16140 

03/14/75 

00% 

8855 

2 

5027 - 16050 

05/16/75 

10% 

5588 


30 -/^ 



ACQUISITION CLOUD 


STATUS 

SCENE ID 
Test Site 3: 

Summer: June - Auq. 

DATE 

COVER 

QUALITY 


1343 - 

16253 

07/01/73 

20 % 

8888 


1361 - 

16252 

07/19/73 

20 % 

8888 


1038 - 

16305 

08/30/72 

20 % 

8888 


1362 - 

16305 

08/30/72 

20 % 

8888 

2. 4 

1703 - 

16175 

06/26/74 


8858 


Fall ; Sept. - Nov. 



1092 

- 16312 

10/23/72 

20% 

8888 


1110 

- 16313 

11/10/72 

00% 

8888 


1452 

- 16291 

10/18/73 

00% 

7828 


Winter: 

Dec. - Feb. 




2, 4 

1146 

- 16314 

12/16/72 

00% 

8888 


1164 

~ 16312 

01/03/73 

10% 

8888 


1182 

- 16313 

01/21/73 

00% 

8888 

2 , 4 

2034 

- 16200 

02/25/75 

00% 

8888 


2016 

-< 16200 

02/07/75 

10% 

5888 

2 

1578 

-16264 

02/21/74 

10% 

8282 


Spring: 

Mar. - May 





1253 

• 16262 

04/02/73 

20% 

8888 


1289 

- 16261 

05/08/73 

00% 

8888 


1236 

- 16320 

03/16/73 

10% 

8888 


1290 

-16315 

05/09/73 

00% 

8888 


1308 

- 16314 

05/27/73 

20% 

8888 

2, 4 

1614 

- 16261 

03/29/74 

10% 

8888 





ACQUISITION 

STATUS 

SCENE ID 

DATE 

CLOUD 

COVER 

QUALITY 

Z, 4 

1614 -^6261 

03/29/74 

10'/ 

8888 


1974 - 16133 

03/24/75 

00% 

8858 

2 

5028 - 16104 

05/17/75 

10% 

8885 

2 

Test Site 4: 

Summer: June - Aug. 
1326 - 16315 

06/14/73 

10% 

8888 

2 

1740 - 16225 

08/02/74 

20% 

8888 


1758 - 16221 

08/20/74 

20% 

8888 

2, 3 

5082 - 16080 

07/10/75 

10% 

8888 


Fall : Sept. - Nov. 


2 

1092 - 16314 

10/23/72 

10% 

8888 


niO - 16320 

11/10/72 

10% 

8888 


1452 ' 16293 

10/18/73 

10% 

8828 


Winter: Dec. - Feb. 




2, 4 

1146 - 16320 

12/16/72 

20% 

8888 


1164 - 16315 

01/03/73 

20% 

8888 

2, 3 

1182 - 16315 

01/21/73 

00% 

8888 


2016 - 16202 

02/07/75 

00% 

5885 

2, 3 

2034 “ 16202 

02/25/75 

00% 

8888 


Spring: Mar, - Ma.y 





1236 - 16323 

03/16/73 

2J% 

8888 


1254 - 16323 

04/03/73 

10% 

8888 


1290 - 16321 

05/09/73 

20% 

8888 



ACQUVSIT70N 

STATUS 

SCENE ID 

DATE 

CLOUD 

COVER 

QUALITY 

2 

1308 - 16320 

05/27/73 

10% 

8888 

2 

1974 - 16135 

03/24/75 

10% 

8858 


5028 - 16111 

05/17/75 

10% 

5588 


Test Site 5: 





Summer: June - Aug. 





1362 - 16315 

07/20/73 

20% 

8888 


1380 - 16314 

08/07/73 

20% 

8888 


1722 - 16235 

07/15/74 

20% 

8888 

2, 3 

1740 - 16231 

08/02/74 

10% 

8888 

2 

1758 - 16223 

08/20/74 

10% 

8888 


Fall: Sept. - Nov. 




2 

1110 - 16322 

11/10/72 

10% 

8888 

2 

1776 16215 

09/07/74 

20% 

5855 


Winter: Dec. - Feb. 




2, 3 

1182 - 16322 

01/21/73 

00% 

8888 


1506 - 16293 

12/11/73 

10% 

8888 

2, 3 

2034 - 16205 

02/25/75 

00% 

8888 


Spring: Mar. - May 

10/18/73 

20% 

8888 


1614 - 16270 

03/29/74 

20% 

8888 


1974 - 16142 

03/24/75 

10% 

8888 


2070 - 16203 

04/02/75 

20% 

8588 

2 

1290-16324 

05/09/73 

20% 

8888 



ACQUISITION STATUS 


1 Imagery on Order 

2 Imagery on Hand 

3 Tapes on Order 

4 Tapes on Hand 



